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An ultrafast paralell data processor is described in which amphtude mode excitations of a charge 
density wave (CDW) are used to encode data on the surface of a l-T TaS2 crystal. The data are 
written, manipulated and read using parallel femtosecond laser pulse beams, and the operation of a 
database search algorithm is demonstrated on a 2-element array. 
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Recently it was shown that the collective amplitude mode (AM) oscillations of a charge-density wave state (CDW) 
can be detected in real time in the quasi one-dimensional semiconductor K0.3M0O3 by femtosecond time-resolved 
pump probe experimentstl. In these experiments, the AM is first excited by an ultrashort laser pump pulse which acts 
as a 5— function-like perturbation, and the resulting reflectivity oscillations due to the modulation of the dielectric 
constant by an AM oscillation are detected by suitably delayed probe pulses. 

In the experiments presented in this letter, we demonstrate the feasibility of using the collective AM excitation of 
a CDW in combination with femtosecond-pulse laser beams to represent data written on the surface of a lT-TaS2 
crystal. We also proceed to perform basic data processing such as storage, manipulation and readout operations. 
Because of the macroscopic nature of the excitation, the operations can be performed at temperatures approaching 
200K with subpicosecond switching and readout times. Finally we demonstrate a parallel search operation using a 
single pulse. 

1 T-TaS2 is a layered nearly two-dimensional metal at room temperature, which undergoes a transition to a CDW 
state driven by a Fermi surface instability at 2kF (kp is the Fermi wavevector) at transition temperatures near 
Tc ~ 200K. The macroscopic CDW state is characterised by two collective excitations: the collective amplitude (AlVl) 
and phason modes (PM), which are described by the solutions of the quantum- mechanical Frohlich Hamiltoniannl. 
Importantly from the point of view of data processing, the AM decay time in this system is quite long, Tp > 10 ps. 
Since this is almost two orders of magnitude longer than the time necessary to excite the mode (~ 100 fs), it 
gives us the opportunity to store data (in the form of phase and amplitude of the AM) and also manipulate the AM 
database in real-time with femtosecond laser pulses. 

In what follows we first demonstrate the manipulation of the AM amplitude and phase using different laser pulse 
sequences, aitpapproach similar to coherent control of coherent optical phonons performed on GaAs and GaAs/AlAs 
superlatticesaQ, Bi filmsO and on crystalline quartzQ. We demonstrate switching of the AM on and off at will, and 
show that we can set and detect the phase and amplitude of the AM to a significant degree of accuracy. In effect, 
we show that the collective amplitude mode excitation can be used as a switchable THz optical modulator whose 
amplitude and phase can be controlled by laser pulses. In the second step we demonstrate the operation of a search 
algorithm with a single-query parallel search on multiple cells defined by multiple laser spots on the sample. 

A schematic experimental setup is shown in Figure 1. A stream of 80 fs pulses from a Ti:Sapphire laser at a 
wavelength of 800 nm is split into a number of beams by Michaelson interferometers in which the time-delay (TD) At 
between the pulses in different arms can be controlled to sub-femtosecond precision. In the standard experimental 
configuration, one or more pulse trains are incident on a spot on the surface of the crystal (a cell), with which the 
amplitude and phase of the AM is written or changed, while a strongly attenuated readout pulse beam is used to 
determine the state of the cells by measuj-ement of the reflectivity modulations as a function of time. The experimental 
details have been described elsewherellI'Ll. 

The reflectivity modulation tp = AR/R at the surface of the crystal caused by a laser pump pulse at t = is shown in 
Figure 2a) (trace A). Typically, the signal can be described by a function of the formtl -iplt, T) = A{t) exp[iw^(r)t-|-(/)] -|- 
i3exp[t/Ts]. The first term describes the AM oscillations, and the second term is the contribution from single-particle 
relaxations. The amplitude of the AM decay is well approximated by an exponential function A{t) = Ao6xp[— i/rp] 
where Tp is the AM decay time. We have found that by tuning the laser wavelength near 800 nm, we can extinguish 
the second term almost entirely, so that only the AM contribution is observed as shown in Figure 2a) . The freauency 
loa{T) is temperature dependent, as shown in Fig. 2b) which agrees well with previous Raman measurementaj. The 
discrete Fourier transform (FT) of in trace A is shown in Fig. 2c). Apart from the AM at uja — 2.42 THz, two 
additional peaks are present, ^2.1 THz and 3.9 THz. These are attributed to coherent phonon oscillations and are 
not of interest to us at presenile. 

In Figure 2a) (trace B) we also show the time-evolution of the AM after the application of a 2-pulse sequence (whose 
timing is indicated by arows). The first pulse sets up the AM oscillatory state with = 0. The second pulse is set to 
have equal amplitude, but is delayed with respect to the first pulse such that its phase is shifted by tt with respect 
to the first pulse, i.e. At — (f) = [2n + ^tt/loa- We see that after the tt pulse in trace B the amplitude of the AM 
oscillations is significantly reduced. From the FT of signal B for i > 2 ps (shown in Figure 2c) it is evident that the 
amplitude of the AM is reduced by more than 90%, while the remaining signal contains significant contributions of 
the coherent phonon oscillations at 2.1 and 3.9 THz, and of course noise. From the point of view of data manipulation 
and computation, we see that the AM has been effectively reset by the tt laser pulse. Considering the fact that no 
spatial or temporal pulse shaping or amplitude matching has been performed in these experiments, the efficiency of 
this process is quite remarkable. 

Similarly, we find that if the second pulse delay is set to be in phase with the first pulse (At ~ (p ^ 2m: /u a), the 
amplitude is approximately doubled at i = by the second pulse, and the total AM amplitude is A — Ai + A2. This 
demonstration shows that A and can be manipulated by the application of an appropriate laser pulse sequence and 
allows any arbitrary superposition of states ip to be represented. 

For the purpose of performing data processing operations, we can divide the surface of the crystal into a set of 
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N = 1...J discrete cells. Since the frequency lua{T) of the AM is temperature dependent (Fig. 2), we can set the 
frequency of different cells with different ujj by adjusting the temperature of each cell. (This can be done with a 
temperature gradient along the x axis for example). This facilitates a tomographic reconstruction of the database 
whereby different cells can be uniquely identified by their frequency, and a Fourier transform (FT) of the readout 
probe pulse which covers the entire array of cells simultaneously gives directly the amplitude (and phase) of each cell. 
In effect it gives a 1-dimensional projection of A and (f> along the x axis. l-l 

To implement a database search, we follow an operation sequence, which is similar to a sequence given by GroverO 
invented for the case of a quantum search. A database search operation is performed in two steps preceeded by a 
setting up operation. First the amplitude Aj and phase of each cell is set up by an appropriately timed " write" 
laser pulse W. This sets up the database ^ such that the information in cell j with frequency uij is represented by 
Aj and . Since the time to perform the entire sequence of operations <C t, we assume that the amplitude is time 
independent Aj{t) ~ Aj and the reflectivity of each cell is thus given by ipj — Aj exp[iujjt + (j)j]. The database is set 
up such that one item is marked by setting its phase to be </) = tt, while the unmarked cells have ^ = 0. 

The next step is to apply 11 pulse on the entire array simultaneously, which amplifies the set item and suppresses 
the rest. In our example this is given by a properly timed laser pulse on the entire array, such that At = (2?! + 1)tt/ujo, 
where loq is the center frequency of the array. The result is that the amplitude of the marked cell is doubled, while 
the rest are suppressed. The reflectivity modulations due to the entire array of cells is then read out by a pulse beam 
_R, which covers the entire array. Thus the data encoded in the surface of the crystal in cells is transferred to the 
reflected readout beam whose intensity is proportional to ^ = Aj expliujjt + (j)j] and contains all the phase and 
amplitude information stored in the array, which can be displayed by a Fourier transform of 

Here we limit ourselves to a two-element array (N = 2), which is sufficient to effectively demonstrate the principle 
of operation. The first cell is heated with a continuous wave argon laser, while the second is not, so that the heated 
and unheated cells have frequencies of wi = 2.38 and lu2 = 2.42 THz respectively. 

A write pulse W first excites the array, with 0i = for cell 1 and 02 = tt for cell 2. The signal ipi in the two 
cells following the operation is shown in Figure 3a). After At = 2mr = 2.09 ps (corresponding to n = 5) we 
apply a 11 pulse. The Fourier components F{uj) at uji and UJ2 before and after the 11 pulse are showa in Fig. 3b). The 
amplitude is obtained by a fit to the data to "i/ ^ + ^2 — Ai cos[[Jit + (j^i] + A2 cos[cLi2i + '^'2]". From the data 
shown in Figures 3a) and b) we see that the amplitude of -01 is approximately doubled, while the amplitude of ■02 is 
extinguished, which, if applied to a larger array, would amount to amplifying the set item and suppressing the rest. 
The identification of the set item thus takes place in a single step. Clearly if > 2, application of the same procedure 
would still yield the same result in a single step, as long as there is only one set item in the database. The same 
search algorithm is also applicable when more than one item is marked, and the amount of information stored in each 
cell is limited by the resolution with which the amplitude and phase can be read out. For a practically useful system, 
to increase N it would be more practical encode the data 2-dimensionally by using position rather than frequency 
to identify each cell, enhancing the inherent parallelism of the processor. The crystal surface is then divided into a 
2D array, where the cells can be constructed of arbitrary size, the lower limit being set by the coherence length of 
the CDW system, ^ 10 nm. Witlt^ome adaptation, the optical elements presently in use for temporal and spatial 
shaping of ultrashort optical pulsealil, can then be used to control the timing of the incident 2D array of pulses. 

To conclude, the demonstrated ultrafast coherent-control data processor and storage device using amplitude mode 
excitations of a CDW is still far from a practical device. However, the parallel processing scheme using a CDW system 
effectively demonstrates a conceptually new approch to ultrafast optical data processing. 

Figure 1 A schematic diagram of the experiments. The pulse timing is controlled by variable time-delay units 
(Michaelson interferometers) TD^. M is a high-frequency acousto-optic modulator. 

Figure 2 (a) The induced reflectivity ip = AR{t)/R of the AM as a function of time t after excitation (trace A). 
The figure shows the raw data. Trace B represents the reflecivity oscillations for a 2-pulse sequence. The timing of 
the second pulse is adjusted so that At — {2n + 1)tt/uja (a tt pulse), (b) The T-dependence of uja- (c) The discrete 
Fourier transforms of ipA and "05 for t > 2.09 ps shown in a). The component at uja is nearly extinguished by the 
application of the tt pulse. 

Figure 3 (a) a real-time display of the parallel search algorithm for a 2-cell array. W sets up the AM oscillation 
with 01 = in cell 1 and 02 = tt in cell 2 at 2.38 and 2.42 THz respectively with equal amplitudes. The operation 
II amplifies the set data item (in cell 2) while suppressing the amplitude in cell 1. The sum amplitude is shown in 
the bottom trace. 

(b) The Fourier amplitudes i^i(w), and phase 0i (expressed both in terms of phase and time-delay) before and after 
the n pulse are shown in the left and right panels respectively. After the 11 pulse, the set item (in cell 2) is amplified, 
while the contents (AM amplitude) of cell 1 are supressed. 
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